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Abstract. Consider a family of bounded domains f2t in the plane (or more 
generally any Euclidean space) that depend analytically on the parameter t, 
and consider the ordinary Neumann Laplacian At on each of them. Then we 

( i) 

can organize all the cigenfunctions into continuous families u\ ' with eigen- 

(?) 

values \f also varying continuously with t, although the relative sizes of the 
eigenvalues will change with t at crossings where A" = \ t . We call these 
families homotopies of eigenfunctions. We study two explicit examples. The 
first example has Qq equal to a square and f2i equal to a circle; in both cases 
the cigenfunctions are known explicitly, so our homotopies connect these two 
explicit families. In the second example we approximate the Sicrpinski carpet 
starting with a square, and we continuously delete subsquares of varying sizes. 
(Data available in full at www.math.cornell.edu/-smh82) 



1. Introduction 

Consider a family of bounded domains fi t in the plane (or more generally Eu- 
clidean space or a Ricmannian manifold) depending on a real parameter t. For each 
domain we consider the standard Laplacian A equipped with Neumann boundary 
conditions. By a homotopy of eigenfunctions we mean a continuous family of eigen- 
functions u t (x) with eigenvalues X(t) (also varying continuously in t e M. 

— Au t — X(t)ut on Q t 
d n u t — on dVt t ^ ^ 

For a reasonable domain f2 t there is a complete orthonormal basis of cigenfunctions. 
We do not expect every cigenfunction to belong to a homotopy. If the eigenvalue is 
simple, then it is easy to see that the (unique up to a constant) cigenfunction does 
belong to a homotopy for at least a small interval in the parameter t. However, if 
the eigenspace has higher multiplicity at a fixed value of t, the best we can expect 
is that there is a special basis of eigenfunctions that belong to a homotopy. 

We face the same situation if we consider a continuous family A t of symmetric 
n x n matrices. The eigenvalues depend continuously on t, but the eigenvectors do 
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not necessarily vary continuously. For example, let 

' t > o 
h 



for Bi, B2 2x2 symmetric matrices with different eigenvectors. One can avoid this 
pathology by assuming that the derivative with respect to t is nonzero, or that the 
dependence on t is analytic. However, if A t — tB where B has distinct eigenvectors, 
then for t = 0, the matrix has eigenspace corresponding to A = containing all 
vectors, whereas only the eigenvectors of B belong to homotopies. 

Observe now that the double eigenvalue of A t at t — is unstable under pertur- 
bation. More specifically, recall that matrices with multiple eigenvalues, considered 
as a subspace of n x n symmetric matrices, have codimension 2. Similarly, those 
matrices with multiple eigenvalues innxn Hcrmitian matrices have codimension 3. 
So, for a generic one-parameter family of self-adjoint matrices, multiple eigenvalues 
are very rare. Moreover, in the self-adjoint case, if two eigenvalues move near each 
other (say, as t increases), the eigenvalues then repel each other. We can see this 
from the following equation (derived below as in the article "When are eigenvalues 
stable?" US]): 

A£ = 4^ + 2£t^¥ (3) 

where u k is the k th eigenfunction (an L 2 normalized column vector), Xj is the j 
eigenvalue, * denotes complex transpose, and ' denotes a derivative in the t variable. 
Note that the eigenvectors and eigenvalues are functions of t. 

To obtain Eq. (|3|, we differentiate the equation A t Uk{t) — Xk(t)uk(t) in the t 
variable to get 

A t u' k + A' t u k = X k u' k + X' k u k (4) 

We then differentiate the identity u* k u k — 1 to get u* k u' k = 0. Then, applying u* k to 
both sides of Eq. Q and using the previous sentence gives 

X' k = u* k A' t u k (5) 

Notice we have also used the fact that 

u* k A t = X k u* k , (6) 

which is just the eigenvalue equation for u k with the adjoint taken on both sides. 
Now, apply u* to both sides of Eq. Q for j =/= k, and apply Eq. (j6| and the 
orthogonality of the eigenvectors to the result to obtain 

u*A\u k + (Xj - X k )u*u' k = 0. (7) 

So, re- writing, summing over j, using the completeness of the eigenvectors (and 
recalling that u*u'j = 0) gives 

<=y:p^ (8) 

Now, differentiating Eq. ^ gives 

x k = u* k A"u k + u* k A' t u' k + (u* k )'A' t u k (9) 
And after applying Eq. @ to Eq. ([9|, we finally get Eq. 
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We now wish to comment further on the significance of Eq. ^ . The repulsion 
of two distinct eigenvalues, or level repulsion, manifests in the second term on the 
right of Eq. pj. That is, if we increase t linearly, and two eigenvalues A& and Afc< 
move close together, the second term of Eq. ([3]) becomes very large (assuming the 
numerator does not vanish). In fact, X'j, has a term of the form l/(Afc — A&'), while 

has a term of the form l/(Afc' — A&). So, the larger of the two eigenvalues is 
forced upward, and the smaller of the two is forced downward. Therefore, Eq. ^ 
accounts for generic eigenvalue repulsion, in this finite dimensional case. 

Now, can these observations carry over to the infinite dimensional setting? We 
are concerned here with eigenvalue repulsion and rarity of multiple eigenvalues. 
Specifically, do these effects occur for compact planar domains which vary with 
a time parameter? To a certain extent, yes. Restricting an infinite dimensional 



operator to finitely many dimensions as in Eq. ( 12 ) below, the repulsion effects from 
the finite dimensional Hermitian case should carry over to self-adjoint operators. 
(For rigorous results on the rarity of multiple eigenvalues, see [Cl ILLI iTe] .) 

Mimicking the finite dimensional case, we reformulate our original problem by 
taking a fixed domain f2 and varying the differential operator. We assume that we 
have a family of diffeomorphisms *5f t : f2 — * £l t that depends analytically on the 
parameter t. Then we can pull back the Euclidean metric on fl t to a Ricmannian 
metric on f2. The pullback of the standard Laplacian on f2 t becomes the Laplace- 
Beltrami operator on J7, which we denote by A t . Similarly, we can pull back the 
energy and measure on f2 t to the corresponding quantities on f2 for the Riemannian 
metric. We note that the Neumann boundary conditions correspond, and also the 
coefficients of A t are analytic functions in t. Let [a, b] be any finite interval in 
[0, oo) such that the endpoints are not Neumann eigenvalues of any A t (to do this 
we may have to restrict the interval in which t varies) . Then the spectral projection 
operator Pt([a, b]) associated with — A t and the interval [a, b] is well-defined and of 
finite rank, and may be expressed as a contour integral in the resolvent 

R(X) = J (X + A^dX (10) 

for a suitable contour T in the complex plane enclosing the interval. The resolvent 
depends analytically on t (this may be seen via the pseudodifferential calculus, for 
example), so the projection operator also depends analytically on t. If the rank of 
the projection operator is N, we can find an iV-dimensional space V of functions 
such that P t ([a,b]) maps V onto its range. Let P t _1 ([a, b]) denote the inverse. Then 
we can transform the eigenvalue equation 

- A t u t = X(t)u t on ft (11) 

into the eigenvalue equation 

- (P t ([a, b])- l A t P t ([a, b])) x t = X(t)x t in V (12) 

via x t = i"t([a, 6]) _1 u t . Now we are in the finite dimensional setting with a family 
of operators depending analytically on the parameter, so N homotopies of eigen- 
functions exist, and these yield N homotopies of eigenfunctions Ut on fl. 

Once we have the existence of homotopies for small intervals of the parameter 
and bounded regions of the spectrum, we can combine them to obtain global homo- 
topies for the whole spectrum. If u to is an eigenfunction corresponding to a simple 
eigenvalue X(t ) of A to , then there is a unique homotopy u t with u t \t=t — u t - If 
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A (to) is an eigenvalue of A to with nontrivial multiplicity, then its cigcnspacc E(to) 
has a basis of eigenvectors which lie in homotopies. In fact one can take this basis 
to be orthonormal. Suppose for example, that A(io) has multiplicity 2, and let 

and u[ 2 ^ denote 2 homotopies with eigenvalues X^(t) and A^ 2 '(i). Then either 
A«(t) = A< 2 )(t) for all * or A«(t) ^ X^(t) for < \t - t Q \ < e. In the first case 
we may take linear combinations of and uf^ to obtain orthogonality. In the 
second case we have the orthogonality of and u-p for < \t — to ! < e (because 
the eigenvalues are distinct), and by continuity it must hold for t — to- 

The deeper question, however, is what is the significance of the homotopies? 
If we have homotopies between domains U and f2i then we have a pairing up of 
the cigenfunctions on Q and f^. What does this tell us about the paired cigen- 
functions? The goal of this paper is to investigate this question experimentally 
We examine two different examples. In the first, something of a warmup exercise, 
we take Slo be a circle and Q\ to be a square. We know exactly what the eigen- 
functions are at both endpoints. For the unit circle, using polar coordinates, the 
eigenspaces typically have multiplicity 2 and are spanned by J m (X m kr) cos m9 and 
J m {X m kr) smm9 where the eigenvalues X m k are the successive zeroes of J' m , indexed 
by k = 1, 2, 3, . . .. Here m = 1, 2, 3, . . .. The exceptional case m — yields a 1- 
dimensional eigenspace spanned by Jo(Aofcr). For the unit square, there is a basis of 
cigenfunctions of the form cosnjxcosTrky for j = 0, 1, 2 . . . and k = 0, 1,2, . . . with 
eigenvalues ir 2 (j 2 + k 2 ). Typical eigenspaces have dimension 1 (j = k) or 2 (j ^ k), 
but higher multiplicity may occur due to coincidences, such as 3 2 + 4 2 = 5 2 + 2 or 
7 2 + l 2 = 5 2 + 5 2 . 

In order to simplify and clarify the discussion we consider the action of the 
dihedral symmetry group D 4 on the eigenspaces. D4 acts on the circle and square 
and all the domains we consider. It is well known that this 8 element group has 5 
irreducible representations, one being 2-dimensional and the others 1-dimensional. 
Each eigenspace either transforms according to one of these representations or splits 
into a direct sum of such pieces, and the homotopies respect the representation 
types. We denote the 1-dimensional representations 1 + +, 1 H — ,1 — h, and 1 — 
— , where the first ± refers to the symmetry or skew-symmetry with respect to 
diagonal reflections, and the second ± refers to the symmetry or skew-symmetry 
with respect to horizontal and vertical reflections. We now describe the eigenspaces 
sorted according to representation type for the circle and the square. 

Circle: 

1 + + : The eigenspace i?+ + (4m, k) is generated by 
J4m{X(4, m )k r ) cos Am6 for m = 0, 1,2, . . . and k = 1, 2,3, . . . 

1 — h : The eigenspace i?i + (4m, k) is generated by 
Jim(X(4, m )k r ) sin4m# for m= 1, 2, 3, . . . and k = 1, 2, 3, . . . 

1 H — : The eigenspace _E^_(4to + 2, k) is generated by 
Jim+2{X(i m+ 2)k r ) cos(4to + 2)9 for m = 0, 1, 2, . . . and k = 1, 2,3, . . . 

1 : The cigcnspacc _E^_(4m + 2, k) is generated by 

Jim+2{X(i m+ 2)k r ) sin(4m + 2)9 for m — 0, 1, 2, . . . and k = 1, 2, 3, . . . 

2 : The eigenspace i?|(2m + 1, k) is generated by 

J2m+i(A( 2m+ i)fcr) cos(2m+l)6> and J 2 m+i(A( 2m+ i)A ; r) sin(2m+l)6> for m = 0, 1, 2, . . . 
and k — 1, 2, 3, . . . 
Square: 
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1 + + : The eigenspace E+ + (2j, 2k) is generated by 
cos n2jx cos 7r2fcy + cos ir2kx cos ir2 jy for < j < k 

1 — h : The eigenspace E^_ + (2j,2k) is generated by 
cos n2jx cos 7r2fcy — cos ir2kx cos ir2jy for < j < k 

1 H — : The eigenspace E+_(2j + 1, 2k + 1) is generated by 

cos 7r(2j + l)x cos ir(2k + l)y + cos ir(2k + l)xcos 7r(2j + l)y for < j < k 

1 : The eigenspace E s __(2j + 1,2k + 1) is generated by 

cos 7r(2j + l)x cos 7r(2fc + l)y — cos 7r(2/c + l)x cos 7r(2j + l)y for < j < k 

2 : The 2-dimcnsional eigenspace E^(2j + 1, 2fc) is generated by 
cos7r(2j + l)xcos7r2fcy and cos7r2fca;cos7r(2j + l)y for < j < k 

In the case of the square, because of coincidences, the full eigenspace associated 
with a given eigenvalue may be a direct sum of some of the spaces described above, 
for example -Ef (5, 0) © £| (3, 4) or i?+_(5, 5) © £+_(7, 1). When we construct ho- 
motopies, we often find that they do not pass through the individual summands. 
This problem does not arise for the circle, since as far as we know there are no 
coincidences among the eigenvalues \ m k- 

In section [2] we describe the homotopies of eigenfunctions for two different ana- 
lytic families of domains joining the circle to a square. In one, the circle expands out 
to a square in horizontal/ vertical orientation, and in the other the circle contracts in 
to a square in diagonal orientation. (The fact that these are not unit squares just re- 
quires a trivial adjustment.) In both cases the homotopies connect the same E s and 
E c eigenspaces. We observe a canonical splitting of two-dimensional eigenspaces at 
the endpoints of the homotopy. We also do not find any two-dimensional eigenspaces 
away from the endpoints of the homotopy. 

In our second example in section [3] we consider a sequence of domains Ai that 
approximate the Sierpinski carpet (SC). Let F±, . . . ,Fg denote the similarities that 
map to square Aq into the squares of the "picture frame" A\ that is Aq with the 
middle ninth removed. Then 

SC = U s i=1 Fi(SC) (13) 

and 

A j = ut^cvo (u) 

In a previous work |BHS] we gave experimental evidence that the spectra of the 
Neumann Laplacians on Aj converge, after the appropriate renormalization, to the 

spectrum of the Neumann Laplacian on SC, with the eigenfunctions also converg- 

(7) 

ing. In other words, there is a complete orthonormal basis u„ of eigenfunctions 

- A W « = (15) 
on Aj with Neumann boundary conditions such that 

lim r j X^ = X„ and (16) 

lim u ( ^\ S c = u n (17) 

j—tOC 

for some r € K. We would then identify u n as the Neumann eigenfunctions of the 
Barlow-Bass Laplacian A on SC [BB . The recently announced proof [BBKTJ of the 
uniqueness of that Laplacian justifies this identification. Combining this uniqueness 
result, the resolvent convergence [B] (see §7), and the relation e tA f(x) = E x [f(X t )] 
between the Laplacian A and Brownian motion X t , suggest a method of rigorous 



proof of Eq. (17 1. The relation e tA f(x) = E x [j(X t )\ is important because it 



6 



STEVEN M. HEILMAN 1 AND ROBERT S. STRICHARTZ 2 



shows that the speeding up of the Brownian motion jB] exactly corresponds to 
multiplication of A by a constant. 

However, even with convergence in Eq. (17 1, there is an obstacle to uncovering 
the dynamics of the eigenvalues and eigenfunctions. It was noted in BHSJ that, 
for a given n, the eigenvalues {An } are not necessarily in increasing order (with 
respect to the original ordering for j = Q). In order to better understand the discrete 
sequence uii of eigenfunctions, we embed them in homotopies. For example, to 
connect the square A with the picture frame A\ , we first remove the center of the 
square. Since points have zero capacity in the square, this does not change the 
Neumann Laplacian. We then gradually enlarge the hole by deleting squares in the 
center with side length going from to 1/3. Similarly, to go from A\ to Ai we 
enlarge 8 holes in the center of the squares F^Aq with side length going from to 
1/9. We carried out the computation of these homotopies up to A 4 . (Note that 
it would also be possible to do a single deformation from Aq to A4 by opening all 
holes simultaneously.) Some eigenfunctions change only slightly during the course 
of the homotopy; we could say that they make minor adjustments to the changing 
domains. However, most eigenfunctions show considerable variation, especially for 
higher eigenvalues. Ideally we would hope to make a reasonable guess from the 
shape of to which u^) it is connected to. We are very far from being able 

to do this. In fact, the following cautionary tale shows some of the obstacles. If j 
and k are divisible by 3, then all the eigenfunctions in any of the E s (j,k) spaces 
are also Neumann eigenfunctions on A\ (when restricted to A\). However, many of 
the homotopies do not connect the restriction to A\ with the original eigenfunction 
on Aq. 



2. Circle to Square Homotopy 

Let L C K 2 denote the line x + y = 1 in the closed first quadrant. In polar 
coordinates L is parameterized by the equation r = 1 /(cos 9 + sin 6) , 9 £ [0,n/2]. 
Then, the map tp(r,6) — (l/(cos 9 + sin 9)r, 9) defined in the closed first quadrant 
takes the arc of the circle of radius 1 centered at the origin to the line L. In fact, 
it also maps the arc of the circle of radius r to the line {(x, y) : x + y = r}. With 
this in mind, let H (t, r, 9) : [0, 1] x R 2 -> R 2 be the homotopy 

H(t, r,d)=( ( (1 - 1) + \ . ) r, 9) (18) 
\ \ cos 9 + sin 9 ) ) 

For 9 £ [0,7r/2], this is the straight-line homotopy between the identity map of R 2 
and ij}(r,0). So, as we vary t, H restricted to the closed first quadrant maps the 
square diamond of side length y/2 to the unit disc. Observe that for t G (0, 1), H 
(with the same restriction) maps the square to a half-ellipse. 

Rotating the domain and range of the map into all four quadrants, we get a 
straight-line homotopy which connects the square diamond to the unit disc and 
which is analytic in the time parameter t on (0,1) (and continuous on [0,1]). If 
we view H(t, r, 9) as an analytic deformation of the standard metric on the unit 
disc, the Neumann Laplacian A t on the disc is a holomorphic family of type (A) 
(see [K]). Then, Theorem 3.9,VII.§4.5 from the same book says the eigenvalues and 
eigenvectors can be represented as functions holomorphic in t. More specifically, 
the eigenvalues arc holomorphic functions and the eigenfunctions can be chosen 
to be holomorphic with a Banach space as their range. In our terminology, we 
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have an (analytic) homotopy of cigcnfunctions. This is the theoretical basis of our 
investigation. 

An analogous process can be achieved for (j>(r, 9) — (v2/(cos 9 + sin 9)r, 9) and 

^^)=((( 1 - f ) + c^f^)^) (19) 

where 4> maps the unit circle in the first quadrant to the square diamond of side 
length 2, and upon reflection as above, F gives a straight line homotopy between the 
identity map on the unit circle and 0(r, 9). Again, we get a straight-line homotopy 
from the unit disc to the square, and with previously mentioned results jK], we 
have a homotopy of eigenfunctions. Note that F and H are not identical. 

Since all of our domains are invariant under the D4 symmetry group, we can 
simplify the eigenfunction computations by reducing to a fundamental domain. 
On this domain we impose appropriate boundary conditions according to the rep- 
resentation type. For the 1-dimcnsional representation, we restrict to the sector 

< 9 < tt/4. Recall that the functions will extend evenly when reflected about 

9 = in the 1 + + and 1 h cases, and oddly in the 1 H and 1 cases. 

When reflecting about 9 — tt/4, the functions will extend evenly in the 1 + + and 

1 H — cases, and oddly in the 1 — h and 1 cases. Note that performing an even 

extension across a ray is equivalent to imposing Neumann boundary conditions on 
that ray. Similarly, the odd extension is equivalent to Dirichlet conditions. 

For the 2-dimensional representation our fundamental domain is the sector 
< 9 < 7r/2, and we impose Neumann boundary conditions on the ray 9 = and 
Dirichlet conditions on the ray 9 = -k/2. Notice that our fundamental domains 
are simply connected. Therefore, as in Theorem 6.4 of |Tej . we should expect no 
two-dimensional cigenspaces along this homotopy (for each individual symmetry 
family) . 

In order for Matlab R2006a to understand our geometry, we must give the 
geometric data to the program in its desired format. To accomplish this task in 
Matlab's PDE Toolbox for a fixed t, we specify polygons and circles. More specif- 
ically, we specify the vertices of the polygon and the center and radius of a circle, 
with standard coordinates in the plane. Once we input this data, we then fix a 
closed interval of the real line in which we look for eigenvalues A. The program 
then uses the Arnoldi Algorithm with function pdeeig to solve the eigenvalue equa- 
tion (— A)u = Xu for A in the specified interval. For this particular homotopy, this 
process yields sixty to seventy eigenvalues within minutes. As our first approxi- 
mation to a smooth homotopy using H or F, we solve for eigenfunctions at times 
t 6 {fc/lOj^o- For now, Matlab's solution method is sufficient. However, other 
solvers based in different programs may be more effective in the future. For in- 
stance, the boundary element method appears much more efficient than the finite 
clement method. 

In the following diagrams, we show the observed eigenvalue dynamics for each of 

the five symmetry families (1 + +, 1 H — , 1 — h, 1 ,2). For convenience, we divide 

all eigenvalues by it 2 . As stated above, the data is collected for t e {fc/10}J£ . As 
t varies, we determine whether or not two eigenvalues cross by inference (see below 
for more details). 

Observe that the eigenvalue plots exhibit eigenvalue repulsion, as in Eq. (§ 
in the finite dimensional setting. When two eigenvalues draw near each other, 
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they "collide" and then veer away rapidly (see Fig. [l9] (a) below). Intuitively, 
when two eigenvalues move near each other, one could imagine a slightly perturbed 
geometry for which these eigenvalues would be the same. Indeed, the eigenfunctions 
also support this intuition. If we examine the eigenfunctions before and after the 
collision, we see that they swap identities. That is, if we have eigenvalues Xj+i > \j 
which collide around t = to, then Uj+i(to — e) looks like Uj(to + e), and Uj(t,Q — e) 
looks like Uj + i(t + e). In the intermediate region (i — e, t + e), we then observe 
two linear combinations of our eigenfunctions Uj+i(tQ — e) and Uj(to — e). This effect 
has been observed before |BSSj and coined the superposition effect. Notice that 
this superposition effect only occurs among eigenfunctions of the same symmetry 
class. 

When a two-dimensional eigenspace occurs (analytically), the superposition 
effect is sometimes not observed numerically. For example, consider eigenfunctions 
iti3 and Uii from the (1++) family on the square. We see from Fig. [ljthat these two 
eigenfunctions separate from the two dimensional eigenspace at t = 1, as t decreases. 
A priori, this two-dimensional eigenspace could split into any two (orthogonal) 
eigenfunctions from a two-parameter family {a\U\3 + (I2U14 : eti, a% G M}. However, 
we see from Figs. [7] and [6] that these eigenfunctions appear to split according to 
the canonical eigenfunctions on the square. This behavior is consistently observed 
for this homotopy. We will see in Section [3] different behavior for that homotopy. 

Additionally, observe the closely bound eigenfunctions in some of our eigen- 
value plots. That is, as t varies in our homotopy, we see two eigenfunctions whose 
eigenvalues remain very close for all t. This effect is puzzling for the authors to 
explain, and it is unclear whether it only occurs in the lower part of the spectrum. 
We observe this effect for eigenfunctions 10 and 11 in Fig. [2] and eigenfunctions 5 
and 6 in Fig. [5] This effect is not observed in our homotopy of Section [3] 

Also note that certain eigenvalues change in similar ways. For instance, eigen- 
values 13,18, and 24 of the circle in the (1 + +) family all have noticeably negative 
slope around t = .5. 

This observation indicates that patterns occur in our eigenfunction homotopies. 
In other words, it seems that sets of similar eigenfunctions behave in similar ways 
in the homotopies. 

In Fig. [7J we depict the homotopy of eigenfunctions of Eigenvalue 13 of the 
circle in the (1 + +) family. In Fig. [8j we show this same homotopy restricted to 
the positive x-axis. 

In Fig[9j we take the same homotopy as in Figs. [7]and[8j but for simplicity we 
only examine the nodal set of the eigenfunction. We approximate this nodal set by 
the inverse image of a small neighborhood of of the form (— e, e). 

In Fig. 10 we depict the homotopy of Eigenvalue 10 on the circle in the (1 ) 

family. In Fig. [TT]we show the homotopy of the corresponding nodal set, as in Fig. 

13 

Observe from Figs. l][5 that the endpoints of the graphs of the eigenvalues give a 
bijection between the eigenfunctions on the circle and those on the square. Since this 
bijection may encapsulate information about our homotopy or the eigenfunctions 
themselves, we examine it. For simplicity, we view the correspondences of the 
nodal sets of the eigenfunctions. Figs. [12|T6| displays pictorially some of these 
correspondences and Table [l] displays all of the data in terms of eigenfunction 



HOMOTOPIES OF EIGENFUNCTIONS AND THE SPECTRUM OF THE LAPLACIAN ON THE SIERPINSKI CARPET! 



labels. For instance, we find nontrivial information about the splitting of the two- 
dimensional eigenspaces on the endpoints (corresponding to t = and t = 1). 

How do we actually determine existence or non-existence of a two-dimensional 
eigenspace as we vary t in an eigenfunction homotopy? Since different symmetry 
families do not interact, it is easy to see when two eigenvalues from different sym- 
metry families cross. For instance, if we combine the (non-interacting) eigenvalues 
of Figs. [T]|5]in a single plot, we will definitely have crossings. But for cigenfunctions 
in the same symmetry family, we cannot rigorously determine from the experimen- 
tal data whether or not two eigenvalues cross, since we are only approximating 
the smooth functions F and G in discrete time. Yet, as we vary the discrete time 
parameter, the graphs of the (lower eigenvalue) cigenfunctions change very slowly. 
So, if two eigenvalues are far apart, and their corresponding cigenfunctions look 
similar as we vary t, then (trivially) these two eigenvalues do not attain the same 
value. For example, in Fig. [l] eigenvalues 1 and 10 on the circle (where t = 0) are 
far apart, and their eigenfunctions look entirely different as we vary t. Therefore, 
these eigenvalues do not cross. 

Subtlety arises when two eigenvalues are very close together. If a multiple 
eigenvalue exists, our numerical methods will not give us two identical eigenvalues. 
This effect is a result of the instability of two-dimensional eigenspaces discussed 
in Section [l] Since we are given approximate solutions to the eigenvalue problem, 
we will get two distinct eigenvalues which are very close together. Measuring the 
difference will tell us the approximate (relative) error of our numerical method. 
Given two nearly identical computed eigenvalues Aj, Aj+i at t = to, we are therefore 
interested in the quantity 



(20) 



For eigenspaces that we know have multiplicity, we find that Ej < 8 x 10~ 5 . Re- 
call from above that, with respect to our eigenvalue data, a collision (see Fig. [l9^a)) 
and a multiplicity are effectively indistinguishable for Ej very small. Therefore, for 
collisions with Ej around 8 x 10~ 5 , we are forced to conclude that a multiple eigen- 
value occurs. Conversely, if £j > 8 x 10 -5 , we know that no multiple eigenvalue 
occurs in a neighborhood of the region we have analyzed. 

In Fig. [T] eigenvalues 12 and 13 on the circle (where t = 0) become very close 
around t = .4 and t = .5. We have observed the discrete homotopy of eigenvalue 12 
in Figs. [7] [8] and [9] In all three figures, we see what appears to be a discontinuity 
in the homotopy around t = .5. 

This apparent discontinuity occurs since the numerical computations we have 
used cannot tell the difference between a two-dimensional eigenspace and an al- 
most two dimensional eigenspace. In other words, for two close eigenvalues, the 
solver combines two distinct eigenfunctions. Our intuition used above can also ap- 
ply here: quite possibly, an eigenvalue collision is just a perturbation (geometric or 
otherwise) of a true two-dimensional eigenspace. Once again, this is a manifestation 
of the superposition effect [BSS . The data shows two linear combinations of the 
true eigenfunctions, instead of a separation into two non-interacting eigenspaces. 
To assuage these concerns, we make higher precision calculations with a higher 
number of time values t around the potential two-dimensional eigenspace. We then 



use the relative error of Eq. ( 20 1 to determine whether or not a two-dimensional 
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eigenspace occurs. If a two-dimensional cigcnspacc is determined to occur, we then 
incorporate the eigenfunction data and the fact that the eigenvalues vary analyt- 
ically (recall Theorem 3. 9, VII. §4. 5 of [K ). If the higher and lower eigenfunctions 
swap appearances around the multiple eigenvalue, then the eigenvalues cross. If no 
such swap occurs, the eigenvalues do not cross. This completes our analysis of the 
spectral dynamics. 

Let us demonstrate the above procedure by example. Along the homotopy 

H for the (1 ) family, we see from Fig. [2] that eigenvalues 10 and 11 on the 

circle (t = 0) get very close together at t = .7. Our depiction of the homotopy of 



Eigenvalue 10 on the circle in Fig. 10 again shows ambiguous behavior at t = .7. 



To make our analysis more precise, we solve the eigenvalue problem for H for 



t G {.6, .62, .64, . . . , .78, .8}. In Fig. 17 we show the tenth eigenvalue for these t 



values, and in Fig. 18 the eleventh eigenvalue. In the time interval [.6, .8], we 



observe that the tenth and eleventh eigenvalues become very close. For instance, at 
t = .72, the normalized eigenvalues are separated by about 1/10, giving a relative 
error E 10 « .00185 > 8 x 10~ 5 . 

So, we conclude that no two-dimensional eigenspace has occurred within this 
time interval. For an example with a two-dimensional cigcnspacc with no apparent 
crossing (i.e. where the ordering of the two eigenvalues does not change), see Figs 
l36l and l37l in Section [3] 



To summarize the different eigenvalue dynamics, we display Figs. 19 and 20 



Fig. 19 shows three different eigenvalue interactions (all from eigenvalue computa- 



tions), and Fig. 20 shows the gaps between each of the respective pairs of eigen- 



values. If we examine the axes of these figures carefully, we see that part (a) has a 
pair of eigenvalues which do not come close enough together for a two-dimensional 



eigenspace to occur. This reaffirms our recent observations of Figs. 17 and 18 
From these same figures we see that the two eigenfunctions swap identities, so with 
our intuition from Eq. ^ in mind, we call this interaction a collision. In part 
(b), we know that a two-dimensional eigenspace occurs at t — .75, but we see from 
Figs. [36| and [37| that the eigenfunctions do not swap identities. Therefore, this is a 
non-crossing. In part (c), we find evidence for a two-dimensional eigenspace, and 
we observe the swapping of eigenfunction identities in Figs. 32 and 33 Therefore, 
this is a crossing. 

We note that the same effect of nearby eigenspaces also occurs on the square. 
Recall that on the square, we have several two dimensional eigenspaces. So, the 
choice of two orthogonal eigenfunctions by our numerical method from such an 
eigenspace amounts to two arbitrary linear combinations of the true eigenfunctions. 
However, as in Fig. [7j we see that, as we approach the square (around t = .9), we 
approach a particular eigenfunction in a two-dimensional eigenspace on the square. 
Since the eigenfunctions vary analytically, the homotopy of eigenfunctions therefore 
selects for us a member of the two dimensional eigenspace on the square. It is 
then interesting to note that, in our observed data, the eigenfunction thus selected 
is always a "natural" eigenfunction on the square, i.e. an eigenfunction of the 
form cos irax cos irby ± cos irbx cos iray, rather than a linear combination of two such 
eigenfunctions. However, in the homotopy of Section [3] we often see that a two- 
dimensional eigenspace on the square splits into two one-dimensional eigenspaces 
in the non- "natural" way. More specifically, the numerical evidence shows a few 



two-dimensional eigenspaces definitely split in a non-canonical way (see Figs. 34 
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and 35 below). For other examples, the evidence is less striking, but non-canonical 
splittings still seem to occur in mostly all examples (bold entries in Table [2| . 

3. Sierpinski Carpet Homotopy 

Let A — {(x,y) S M 2 : max{|cc| , \y\} < |} be the unit square centered at the 
origin. Let {/j}f =1 denote the usual eight afhne linear maps used in the construction 
of SC, and let 

8 

A i =\Jf i A j - 1 , (21) 
i=l 

for j > 0, j € Z. Recall that the Sierpinski carpet SC is defined (as in Section [T]) 
by 

OO 

SC=f] A l . (22) 

i=0 

We would like to deform A\ continuously into Aq, and then iterate this procedure 
for higher Aj. (The choice to make a map Ai — ► Ao instead of Ao — ► A\ is for 
convenience. Also, as stated in Section [T] note that the removal of a single point 
does not change the Neumann Laplacian.) To construct these deformations, we 
follow the procedure of Section [2] Let 

C= {(x,y) GR 2 : x>0, J > tan" 1 (y/x) > -^}u{(0,0)} (23) 

so that C is the right "quarter" of the plane. Then let L = A\ H G be the right 
"quarter" of the picture frame A\. Now, define G(t,x,y): [0,1] x L ^ A nC by 

G(t,x,y)=(^tx-^t + lj-(x,y). (24) 

Observe that G is a straight line homotopy from the identity map to the map 
(x,y) h-> ((3/2)x - (1/2)) ■ (x,y). As in Section [2] if we rotate the domain and 
range of the map three times by increments of tt/2 radians, we obtain a continuous 



map G : [0, 1] X A% — * Aq. This map is depicted in Fig. 26 To extend G to the 



remaining Aj, j > 0, define Gj(t,x,y): [0, 1] x Aj + i — > Aj inductively by 

8 

Gjit, x,y) = U fi o Gj-i(t, fr\ x , y)). (25) 



i=l 



In words, we apply Gj-i to each of the eight images {/iAj_i}f =1 used to define Aj. 
So, for a fixed j, we use fi : Aj-\ — > Aj, so that /r : fiAj-x — * Aj-i is the (well- 



defined) map in Eq. ( 25 1 . Finally, since Go = id on <9Ao C Ai for all t, the map Gj 
is well defined (a priori, each Gj(t, •, •) has two definitions for certain (x,y) € -Bj). 
The maps G\ through G3 are depicted in Figs. [27| through [29| respectively. 

As in Section [2] we see that, for each j, the eigenvalues and eigenvectors of Gj 
can be represented as functions holomorphic in t. However, since the domains of Gj 
and Gj+i only meet continuously, the eigenvalues only vary continuously around 
im(Gj(l,-,-)) = im (G J+ i(0, ■, ■)). We further note that Gj{X,x,y) maps Bj to 
Aj+i as expected. Finally, to stay consistent with Section [2] we refer to Gj(t,x,y) 
below as Gj(l — t, x, y) as defined in the previous few paragraphs (we have already 



used this convention in the caption of Figs. 19 and 20) 
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As above, our domains {-Bjj^Lo are invariant under the -D4 symmetry group. 
In our computations, we restrict each Bj to the sector 0<#<7r/4in polar coor- 
dinates, for the 1-dimensional representations. We impose our four combinations 
of boundary conditions on the intersection of our domain with the rays 9 = and 
9 = it/ 4. For the 2-dimensional representation, our fundamental domain is Bj n L. 
We impose Neumann boundary conditions on the ray 9 = 7r/4 and Dirichlet condi- 
tions on the ray 9 = — tt/4. 

For all fundamental domains, we place Neumann boundary conditions on all 
remaining edges. Notice that once again, our fundamental domains are simply 
connected for j = 0. Therefore, as in Theorem 6.4 of |Te| . we should expect no 
eigenvalue crossings along this homotopy (for each individual symmetry family). 
Since eigenvalue crossings do in fact occur, this suggests that there is something 
exceptional (i.e. non-generic) about our domains. The appearance of these crossings 
could be related to the quantum-classical correspondence principle. More specifi- 
cally, since the billiard flow is not ergodic in phase space, this may suggest a dense 
number of crossings. 

In Figs. [21] to [25] we display the observed crossing of eigenvalues for the five 

symmetry families (1 + +, 1 H — , 1 — h, 1 ,2) with eigenvalues divided by it 2 . 

We collect the data of G 3 {t, •, •) for j = 0, 1,2,3 and for t € {k/10}l° =0 - Observe 
that, in the bottom of the spectrum , the incidence of crossings decreases as j 
increases. More specifically, the incidence of non-crossings grows exponentially. 
Also observe that, indeed, the eigenvalues of Gj and Gj+i meet in a continuous but 
non-differentiable manner. 

In Figs. [26]|29l we show the homotopy of eigenfunctions of Eigenvalue 4 in 
the (1 + +) family for Go through G3. This eigenfunction exhibits the typical, 
expected behavior. The function changes a great deal for Gj with j small, and 
then as j grows large the eigenfunction maintains the same appearance. However, 
the eigenvalue continues changing, as expected. 

In Fig. [30] we show the nodal homotopy corresponding to Fig. [26] With a low 
eigenvalue, this eigenfunction exhibits minimal irregularity during the homotopy. 



In Fig. 31 we show the homotopy of 26 restricted to the line {(x, y) £ K 2 : < x < 
1/2}. _ 

It is interesting to note how the two dimensional eigenspaces split for Go(i, •, •) 
as we move t from to some small e > 0. In Section [2] we saw a "canonical 
splitting" for the given homotopy. However, the eigenfunctions for Go(e, •, •) for two 
dimensional eigenspaces do not resemble the natural eigenfunctions of Gq(0,-,-). 



For example, in Figs. 32 and |33[ we see the two-dimensional eigenspace on the 



square split with an apparent discontinuity. In Figs. 34 and 35 we see a (non- 
canonical) splitting of the two-dimensional eigenspace on the square. If this were a 
canonical splitting, we would have the (10,0) and (8,6) eigenfunctions. However, 
there is no direct rcscmblcnce of our computed eigenfunctions and the canonical 
ones. In fact, around t = A, these eigenfunctions cross again, since the relative 
error between the eigenvalues is about 6.4 x 10~ 5 < 8 x 10~ 5 (recall the discussion 



around Eq. (20 1. Finally, in Figs. 36 and 37 we show the homotopy around an 



(analytically determined) two-dimensional eigenspace at t = .75. It is not hard to 
find infinite sets of such eigenspaces for any rational t € (0, 1). 



homotopies of eigenfunctions and the spectrum of the laplacian on the sierpinski carp eos 

4. Conclusions 

We now summarize our results. In both of our homotopies, we preserve the 
D4 symmetry of the domains. This allows us to split the cigcnfunctions on these 
domains into five non-interacting families. That is, the behavior of cigcnfunctions 
in one symmetry family along the homotopy does not affect the other families. The 
following is then clear: eigenvalues of one symmetry family definitely cross with 
eigenvalues of other symmetry families. However, the numerical evidence leads us 
to speculate on the individual symmetry families as follows. When we deform the 
circle to the square, we observe no intermediate two-dimensional eigenspaces. How- 
ever, when we deform the square to the picture frame (a square with a smaller square 
removed from the center) we do find intermediate two-dimensional eigenspaces. In 
fact, analytically one can easily find an arbitrarily large number nontrivial multi- 
plicity eigenspaces for a fixed rational time. 

When we observe our domains, we see that the deformation of the circle to the 
square is a convex diamond-shaped region for all t <G [0, 1]. In fact, for t £ (0,1), 
the region is a diamond with four edges, each of which is a half-ellipse. However, 
as we map the square A to the picture frame A\ via H , we always have a rational 
polygon (a polygon with angles which are rational multiples of 7r). Mapping Ai to 
Ai + i via Gi always yields a rational domain with a finite number of flat edges (which 
is now multiply connected for i > 1). Recall: (a) after reducing multiplicities due 
to symmetry, we observe eigenvalues of nontrivial multiplicity for Go but not H. 
Also, (b) for Gi we observe non-canonical splitting of two-dimensional eigenspaces 
of the square, and for H we only observe canonical splitting. 

With respect to observation (a), it is tempting to make a connection between 
the appearance of eigenvalue multiplicity and ergodicity of the billiard flow, but 
such a connection is unclear. Such considerations seem related to the quantum- 
classical correspondence, and some conjectures discussed in Sa . In particular, the 
author mentions that, for the bottom part of the spectrum, the eigenvalues for 
integrable billiards may follow a Poisson distribution, while the eigenvalues for a 
chaotic system may follow a GOE (Gaussian Orthogonal Ensemble) distribution. 
With respect to observation (b), the authors have no clear explanation. 

We have just discussed possible connections between ergodicity of billiard flow 
and eigenvalue properties. Such connections are collected under the so-called 
quantum-classical correspondence. Perhaps extending these connections further, 
we discuss the localization and near-localization (or scarring) of cigcnfunctions. 
For the standard Laplacian on the Sierpinski Gasket (SG), many cigcnfunctions 
are localized [St]. Indeed, it appears that for many symmetric PCF fractals, the 
symmetry of the fractal and the capability to break it into disjoint pieces by re- 
moving finitely many points contributes to the localization of eigenfunctions. In 
contrast, we have SC, which is a symmetric non-PCF fractal (intuitively, this means 
it cannot be broken up into smaller pieces by removing finitely many points). 

We now ask: do the "stronger" connectivity properties of SC forbid localized 
eigenfunctions? We believe the numerical evidence from this paper suggests, yes. 
Indeed, no localized eigenfunctions appear in the computations. If the outer ap- 
proximation method of |BHS] actually works, then this observation finds strong 
support. For possible counter-evidence, let us consider the results of |HH| . which 
imply that any rectangular portion of a domain almost always exhibits eigenfunction 
scarring (that is, having mass concentrated in the domain in a non-equidistributed 
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way) . This would imply that the number of scarred modes for our domains should 
increase exponentially as j — > oo as in Eq. (17 1. However, it is unclear whether a 
scarred eigenfunction for a fixed j would remain scarred as j increases. If the outer 
approximation method of [BHS holds, then the eigenfunction scarring may have 
repercussions for the eigenfunctions on SC , even though strict localization may not 
occur. 

Although it may seem unfortunate that we are not able to arrive at any coherent 
description, even just conjectural, relating the behavior of eigenfunctions joined by a 
homotopy, there is also a more optimistic conclusion: the connection is rather com- 
plicated, and requires further investigation. We hope that the data we present will 
inspire further investigation. (Data available in full at www. math. Cornell . edu/~smh82) 
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Eigenvalue Correspondences 

Eigenvalue (1 + +) Family (1 ) Family (1 - +) Family (1 + -) Family (2) Family 

Number on 
Circle 





Circle 


Square 


Circle 


Square 


Circle 


Square 


Circle 


Square 


Circle 


Square 


1 


( 1, 0) 


( 2, 0) 


(0, 4) 


( 3, 1) 


( 0, 2) 


( 1, 1) 


( 0, 2) 


( 2, 0) 


( 0, 1) 


( i, o) 


2 


( 0, 4) 


( 2, 2) 


( 1, 4) 


( 5, 1) 


( 1, 2) 


( 3, 1) 


( 1, 2) 


( 4, 0) 


( 0, 3) 


( 2, 1) 


3 


( 2 , 0) 


( 4, 0) 


(0, 8) 


( 5, 3) 


( 0, 6) 


( 3, 3) 


( 0, 6) 


( 4, 2) 


( 1, 1) 


( 3, 0) 


4 


( 1,4) 


( 4, 2) 


( 2, 4) 


( 7, 1) 


( 2, 2) 


( 5, 1) 


( 2, 2) 


( 6, 0) 


( 0, 5) 


( 3, 2) 


5 


( 0, 8) 


( 4, 4) 


( 0, 12) 


( 7, 3) 


( 1, 6) 


( 5, 3) 


( 1, 6) 


( 6, 2) 


( 1, 3) 


( 4, 1) 


6 


( 3, 0) 


( 6, 0) 


( 1, 8) 


( 7, 5) 


( 0, 10) 


( 5 , 5 ) 


( 0, 10) 


( 6, 4) 


( 2, 1) 


( 4, 3) 


7 


( 2 , 4) 


( 6, 2) 


(3, 4) 


( 9, 1) 


( 3, 2) 


( 7, 1) 


( 3, 2) 


( 8, 0) 


( 0, 7) 


( 5 , o) 


o 
o 




y o, 4j 


( A °) 


( Q 0\ 
\ **J "/ 








K °, ) 


y l, o) 




9 


( 0, 12) 


( 8, 0) 


( 0, 16) 


( 9, 5) 


( 0, 14) 


(7, 5) 


( 0, 14) 


( 8, 4) 


( 0, 9) 


(6, 1) 


10 


(1,8) 


( 8, 2) 


( 1, 12) 


( 11, 1) 


( 4, 2) 


(9, 1) 


(4, 2) 


( 8, 6) 


( 2, 3) 


( 5, 4) 


11 


(3, 4) 


( 6, 6) 


(4, 4) 


( 9, 7) 


( 1, 10) 


( 9, 3) 


( 1, 10) 


( 10, 0) 


( 3, 1) 


( 6, 3) 


12 


( 5, 0) 


( 8, 4) 


(3, 8) 


( 11, 3) 


( 3, 6) 


( 7, 7) 


(3, 6) 


( 10, 2) 


( o, ii) 


( 7, 0) 


13 


( 2, 8) 


( 10, 0) 


( 0, 20) 


( 11, 5) 


( 5, 2) 


( 9, 5) 


(5, 2) 


( 10, 4) 


( 1, 7) 


( 7, 2) 


14 


( 0, 16) 


( 8, 6) 


(5, 4) 


( 13, 1) 


( 0, 18) 


( 11, 1) 


( 0, 18) 


( 10, 6) 


( 2, 5) 


( 6, 5) 


15 


( 1, 12) 


( 10, 2) 


( 2, 12) 


( n, 7 ) 


( 2, 10) 


( n, 3) 


( 2, 10) 


( 12, 0) 


( 3, 3) 


( 7, 4) 


16 


( 4 , 4 ) 


( 10, 4) 


( 1, 16) 


( 13, 3) 


( 1, 14) 


( 9 , 7) 


( 1, 14) 


( 12, 2) 


( 4, 1) 


(8, 1) 


17 


( 6- 0) 


( 8, 8) 


(4, 8) 


( 13, 5) 


( 4, 6) 


( 11, 5) 


(4, 6) 


( 12, 4) 


( 0, 13) 


( 8, 3) 


18 


(3, 8) 


( 10, 6) 


(6, 4) 


( 11, 9) 


( 6, 2) 


( 9, 9) 


(6, 2) 


( 10, 8) 


( 1, 9) 


( 9, 0) 


19 


( 0, 20) 


( 12, 0) 


( 3, 12) 


( 13, 7) 


( 3, 10) 


( 13, 1) 


( 3, 10) 


( 12, 6) 


( 2, 7) 


( 9 , 2 ) 


20 


(5, 4) 


( 12, 2) 


( 0, 24) 


( 15, 1) 


( 0, 22) 


( 11, 7) 


( 0, 22) 


( 14, 0) 


( 0. 15) 


( 7, 6) 


21 


( 2, 12) 


( 12, 4) 






( 2, 14) 


( 13, 3) 


( 2, 14) 


( 14, 2) 


(3, 5) 


(8, 5) 


22 


( 7. 0) 


( 10, 8) 






( 5, 6) 


( 13, 5) 


(5, 6) 


( 12, 8) 


(1,11) 


( 9, 4) 


23 


( 1, 16) 


( 12, 6) 


















24 


( 4, 8) 


( 14, 0) 


















25 


( 6, 4) 





















Table 1. Table of Eigenfunction Correspondences: For the 
Circle, (a, b) <-> Jh(A„ ja r) cos(27r665). For the Square, (a, 6) <-» 
cos(7raa:) cos(7r6y) ± cos(7ray) cos(7r6a;), (+ in columns 1,3, — in 
columns 2,4). Two dimensional eigenspaces in bold. 
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Eigenvalue Correspondences 



(1 + +) Family (1 - -) Family (I + -) Family (1 - +) Family (2) Family 



Eigfcn 
Square 




<4i 


A 2 


^3 


4 j 


Eigfcn 
Square 


A 




A? 


A3 


Aa 


Eigfcn 
Square 


^0 


A i 


Ai 


Az 


Ai 


Eigfcn 
Square 


Ao 


^1 


A2 


^3 


Ai 


Eigfc 11 
Square 


A 


^li 


A2 


^3 


Ai 


( 2, 0) 


— 




4 






11 11 


4 


4 


4 


4 


4 


( 2, 0) 


4 


4 


4 


— 


— Y 


M 11 


— 


4 


I 


1 


— j- 




1 


1 


4 


4 


4 


(2 2) 


2 


2 


2 


2 


2 


( ^ 11 


2 


2 


2 


2 


2 


( 4 0) 


2 


2 


2 


2 




/ '< 11 


2 


2 


2 


2 




< 9 1 1 


2 


2 


2 


2 


2 


f 4 , 0) 












I i>, 3J 












1 A 01 












( 3, 3) 










,^ 














( 4 21 


4 


4 


4 


4 


4 


(7 11 


4 


4 


4 


4 


4 


( 6, 0) 


4 


4 


4 




4 


1 5 11 


4 


4 


4 


4 


4 


{ 'i 21 


4 


4 


4 


4 


4 














( 7 11 












( 6 2) 












( ^ 11 












1 A 11 












(6 0) 


g 


g 






5 


( 7 ^1 


„ 




g 


g 


g 


( 6 4) 




6 


_ 


6 


6 


t 7 11 


g 


g 


g 


g 


6 


( A t\ 


6 


7 


7 


7 


7 


( 6, 2) 






7 




7 


1 q 11 


7 


7 




7 


7 


( 8, 0) 


7 


7 


7 


7 


7 


( 5, 5) 


7 


7 


7 


7 


7 


( 5> 0) 


7 






g 






g 


g 








1 q 11 


g 


g 


g 






(8 2) 












( 7, 3) 












( 5 2) 




g 


g 






( i 0) 


9 


9 


11 


11 


11 


(9, 5) 


9 


9 


9 


9 


9 


( s! 4) 


9 


9 


9 


9 


9 


(7, 5) 


9 


9 


9 


9 


9 


(6, 1) 


9 


9 


9 


9 


9 


( 8, 2) 


10 


10 


9 


9 


9 


( n, 1) 


10 


10 


10 


10 


10 


( 10, 0) 


10 


10 


10 


10 


10 


(9, 1) 


10 


10 


10 


10 


10 


( 5, 4) 


10 


10 


11 


11 


11 


( 6, 6) 


11 


11 


10 


10 


10 


(9, 7) 


11 


11 


11 


11 


11 


( 8 > 6 ) 


11 


11 


11 


11 


11 


(9, 3) 


11 


11 


11 


11 


11 


( 6, 3) 


11 


11 


10 


10 


10 


( 8, 4) 


13 


13 


13 


13 


13 


( 11, 3) 


12 


12 


12 


12 


12 


( 10, 2) 


12 


12 


12 


13 


13 


(7, 7) 


12 


12 


12 


12 


12 


(7, 0) 


12 


12 


12 


12 


12 


( ') 


13 


14 


14 


14 


14 


( 11- 5) 


13 


13 


13 


13 


13 


( 10, 4) 


13 


13 


14 


14 


14 


( 9, 5) 


13 


13 


13 


13 


13 


(7, 2) 


13 


13 


13 


13 


13 


( ') 


14 


13 


13 


13 


13 


( 13, 1) 


14 


15 


15 


15 


15 


( 10, 6) 


14 


14 


13 


12 


12 


( 11, 1) 


14 


14 


14 


14 


11 


( 6, 5) 


14 


14 


14 


14 


14 


( 10, 2) 


15 


15 


15 


15 


15 


( 11, 7) 


15 


14 


14 


14 


14 


( 12, 0) 


15 


15 


15 


15 


15 


( H, 3) 


15 


15 


15 


15 


15 


(7, 4) 


15 


15 


15 


15 


15 


( 10, 4) 


16 


16 


18 


19 


19 


( 13, 3) 


16 


16 


16 


16 


16 


( 12, 2) 


16 


16 


16 


16 


16 


( 9, 7) 


16 


16 


16 


16 


16 
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Table 2. Tabic of Eigcnfunction Correspondences: For the 



Square, (a, 6) <-> cos(7ra;r) cos(nby) ± cos(7ray) cos(7r6a:), (+ in 
columns 1,4, — in columns 2,3). Two-dimensional cigenspaces 
in bold. Question marks denote blatantly non-canonical splitting 
of a two-dimensional eigenspace 
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Figure 1. (1 + +) Eigenvalue Summary Plot. Circle (t = 0), 
Square (t = 1) 




FIGURE 2. (1 H — ) Eigenvalue Summary Plot 
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FIGURE 3. (1 - +) Eigenvalue Summary Plot 




FIGURE 4. (1 ) Eigenvalue Summary Plot 
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Figure 5. (2) Eigenvalue Summary Plot 
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l.i . N(HmVa): ,50,0430 1-0.3. NotmVsl: ,+6,3905 1-0.9, NcnmVal: ,46.5374 



♦ 







1=0 7, NwmVt* ,44.7460 1-0.6. NwmVal: 4£ 1-0.5, NcrmVal ,41 .0326 



MU, NonWal ,39.2093 1-0.3. NotmVal: ,37.4615 1-0.2, NcnmVal: .35.S374 

4fe Jfc 



1-0.1. NafmVEJ: .34.3S13 t-O, NormVal: ,33.0307 



Figure 6. (1 + +) EigVal 14 Homotopy 
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l.i . NwrnVrf: .50 0434- 




f-0.7, NwmVd: 41.1862 
V0.4, NotmVul: ,33.3945 

# 

1-0.1 , NocmV^: .33.2235 




t-O.S. NotmVal: 47.103S 




1.0.6 NwmVal: .36 3*07 




1.0.3 NmiiiVjI: ,34^359 




Wt, hormVal: ,32.0286 




1=6.3, NormVal ,44.033* 

# 

1-0.5, NcrmVal ,36.2063 




IdO.£, NcnmVal ,34.1503 




Figure 7. (1 + +) EigVal 13 Homotopy 
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1=1 NcimVal^, 0294 



0.0, NonmVal 3S 5354 




WQJB, N«mVnl:37,504 




t-0.7, MafmV^:3S.S061 



1=0.6. NcrmVaJ3SJ56 



1-0.5. MofmVd35.E27S 




t-OA, N«mVal:33.99£2 



1-0.3, HormVd:3£.Q435 



t-Oi, Not mVal:3Q.3(}9 




t-0.1,NarmVd:E8.7818 



1*0. NormVal:27-5029 




Figure 8. (1 + +) EigVal 13 Homotopy, Restricted to a Line 
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1-1 , HormVal: ,4O.0M4 1*0.8, NormVal: 33.5554 1*0.3, NormVal: ,37.594 




1*0.7, NomVal : ,365061 1-0 .6, NormVal: ,35.35 2 1=0.5, NormVal: ,555276 




1*0.4, NamWal: ,335912 1*0.3, NormVal: ,32.0*65 1*0.2. NormVal: ,30.309 




Figure 9. (1 + +) EigVal 13 Nodal Homotopy 
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l.i . NormVal: &i 0653 1-0.3 . NormV*| : 5$ 1-0 .9, NcnmVal: ,55.0372 



♦ 






M)7, NormVak,S3.0i£4 1-0.6. No«mV*l: ,4^C0S3 1-0 .5, NcrmVal ,4$.e57& 



M) 4, Nor mVi* ,43.9053 1-0.3 Normal: 41 4205 M) .2, NcnmVal: ,39.206* 

jfe ^ftt 



1-0.1 . NofmVEJ: .37.3674 t-0, MormVal: ,35.6243 








Figure 10. (1 ) EigVal 10 Homotopy 



HOMOTOPIES OF EIGENFUNCTIONS AND THE SPECTRUM OF THE LAPLACIAN ON THE SIERPINSKI CARPF2F 



c^*p4^jif*« do. i *^ *v ijp r 



1-1, NcnmVak ,01.0685 1*0.8, Norm VaJ: ,53.7969 NonmVal ,56.C37£ 




1*0.7, NomVal: ,53.01 24 1*0.6, McrmVaJ: ,49.6933 1-0.5, MormVal: ,46.6576 




1*0.4, NamWal : ,435053 1*0.3, MormVal : ,4 1 .4205 1-0.2, MormVal: .39.2064 




1*0.1, NomVal: ,37.2674 1*0, NcrmVal: ,35.6243 




Figure 11. (1 — 



-) EigVal 10 Nodal Homotopy 
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V^SlNcrmval:1.43E V^SlNormval:2 

o <§> 

V8W!2Ncnnval:£.e6« VsWSSN&rmvaM 

VilfSNcnn valvar Vaks3Nonmval:S ODi 

® # 

V^*4Ncrmval:3.73£ Val#4Normval:1 
ft] 



V«l*SNemiv*l:9.4M VaWSNcnwaM & 



VJ«6Ncrmval: 10.43 ValftBMormval: 1 3.01 



VsW7N*rmval 1«3 Vala7Nor™val:20.Cii 
v; *::Mmmu- il:1T.« Val#8Noffmral:2&.0i 
ValS9Normval:19.5£ ValS@rSomwal:32.£}£ 



VaJSIlNafmvaJiSjEJ Va]#1 lNcrTnva]:38.02 
V3*Sl2Neim«5J27.S ValSI SHomwaliJO 03 

HI 



ValS14Nofmval:33.0E Val#14Ncrmval:50.05 



ValsiSN&fmval 35 63 Val#15Ncrmval:S2.CM- 



ValS16Nofmval:37^7 Vd#1SNcrmval:5£.06 



Valai 7No«mval 39 02 Vd#17N«mvdj64.07 



VsW13f^4tmv«J45.7 Val#1EMcrmval:Se.09 



VaIS1SNarmvaI:50.07 ValSI 9N<hiihh^:72.09 



VsHGWlixm** 50,9 Val#MNflnnv*l:74,M 



Figure 12. (1 + +) Eigenfunction correspondences 
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VaifflNcrmvaJ:2.£65 Valfi1Normval:5 

VaW2N*imval:S.ro2 V*WI2Nflrmval:1S 
Vi**3Ncrmval:9.452 VaW3Normval:17 



Val#4M(Hmvd:18.3 Val*f4NDrmval:25.Q1 



VnW5Ncrmval:19.S2 V*W5Normv»l:29.01 



VaW6Naimvd:20.2 Val*6NDrnwal:37.QE 



VaH7Ncrmval:£3.S3 Vals"Normv*bU.M 



Vt**6Ncmival:3fi.W V*l*eNorwal:4i.*4- 



VJff9Norrnval:33.0B ValS9 Normal: 53.05 



ValSlON(Kmv;il 3S 6 H V:<H!1 ONurmual 6 1 . 7 




Val*11 Norm val :37.3 a ValaUNomwali&S 07 



Val*12Nsrmvali4S.7 W*f2Num V fl:6S.07 



VMl«nN«,in»a):50.07 Val*T3Nom*8l 73.09 



Valid 4Monnval :50. 9 VaHf 1 *N 'acmval£ 5. 1 2 



V ; .la l SNemvaJiSI .93 Val*t 5No«n»«J 35.13 



Valid 6Normval:54.39 Val#16Normv.al:3S_15 



Val*1 7INoimva):61 .33 Val# 1 7N araival37. 1 7 



Vdl* l 8M<ji rfi ;$6 ,44. Val#1 aiMarra va!:1 01 .2 




Valid aNormval:Ga.ag Val ttf 9Narmval:1 09.2 



V9l*20NoiTnv8l:70A» Val*2QNcrmval:113.2 



Figure 13. (1 ) Eigenfunction correspondences 
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ValSlMomival:0JS452 Val#1 Normvd 1 

X 

Valfl2N4irnval^.SS7 V«|4SNoimv&l:5 
VnW3Nflnnval:5.7(£ VaUI3Nofmval:9.001 



ValS4Mcrmval: 10.07 VaBMNaf mv^:1 3 

4*? 



VaW8Nflrmval:13 96 VaWSNormval 17 01 



VaUBNcrmval: 14.04 W-HVn ~nn C ' 



VtW7Nflrmval:17,5S Va1ff7hormval:25.01 



Vl*»Nflmival:S3 63 ValSeNormval:£9.a£ 



VaiffSMcrm val : 25.37 Val« Moi - val:37.0G 



VaWH lNomiYal:2742 VeJSI 1 Norni val:45 04 



VaU1SNennval:3S.» VaUISNoimvalnS 04- 

I 

VaU 1 3Normval:3».e 1 Val#1 SNorm val :S3.05 



Val#14toomwail:41.15 Valsi4No«ni'al 61 ,0& 



Valsi 6Nof m val -4+ 7S Val* 1 eNcrm val :63.tW 




ValS17Narmval:4E.7e V'alS17Natmval:73.1 



Valff18Narmval:B£.14 VaUM3U*mval 31 IS 



ValS1SNwmval:57.2T ValS13HasmwaJ*5.13 



Val*10NwmvalS7 0fl Val#10Normval:41.aa| ValffEONafmvaliSSjES Val*aoNormval:35.13 



Figure 14. (1 — h) Eigenfunction correspondences 
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Va if-1 Normwal:0.9452 Val H 1 Noting :£ 



Val«SNoaiwtiJ.4. Si7 VaW2Natmwal:8.M1 







Val#3NM™v*;3.7C2 Va S2\cmval:10 



Val#4N«mval:10.07 VtWtNcrfnwal 1 S 1 



VaKfSNomxgJ 13.3& VaH5Normval:2D.D1 



Val*6H«mval 14,04 Vala«Nor™val:2&,9i 




Val#7NaniiYal:17.53 Vdff7f*)n™al:32.02 



VaI#3N«mval£3.63 Val«BHomwal:34.(ie 



Va # 1 C \ di -f wal:27.09 Val#1 ONorm val :50.04 



ValS1 1 Norm val:27 J2 VaW 1 1 Ncrm val :SC .04 




VaWi2Nwmval:35,2£ Val#1£Normval:5£.Q4 



Val*13Nonn™i:38.ei Val#13NorTnval:58.W 



ValSUNorm val:41 .1 5 Val Jl 1 4Ncrm val :6S,tW 



V*Wl3No*mv*l 41 47 Val#15NorTnval:72.09 



ValSI eNo* m val 44 73 Valfll 6 \ c — val:74.09 




ValS17Normval:48.7S VaWI 7Normval:3C.1 1 



Valff18Normval:5£.14 Val*1«Nflrmval:32.1 1 



VaW1SNwmval:57 2T VaUMSNcrmvakSO 14 



Valff20Nofmval:5S J e5 Vsf&tWewiK'd 



Figure 15. (1 H — ) Eigenfunction correspondences 
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VsMIINfXITwal 0.3455 ValilMcmiYal:0_5 



Val#2ti«mval:1.783 Valfl2Normv«l:2.8 



Va\S3tkrmv£:Z£S Vals3Normvd:4jS 
Val#INo™wills4J71 VaW4NatrTwal:6.501 



7d#5>&nnvaJ:6i1 V^SN^mval 9,502 



Viil*e\i:':~vfi VeiISG\c — va "2.5 



Val#7M(Hmval:7.45S VtH7Ncrm val : 1 2 5 




Val#3Non™ai:1 1 22 V a WSNe«n*a! M.3 



Va|i9fJe*™*a* 1 1 .03 ValffSNormraJ: 1 3.51 
Vafcl ONcnnval: 13.05 WslONoflnvg 20.31 



ValSI 1 Nairn val: 13,69 VaWI lNcrmval:22 51 



4§k 



Valfll2N<Hmval 10.07 Va1#12Ncrmva1:St.51 

- 



ValSlSNofm val 10 93 ValSI 3Nonnual:26.52 



4- 



ValS14Narmval:1Sj8S ValSlJNoffnvgl 30.32 
VaUn3Ne*mval2i 30 Va|4MSNflnnval:32,52 



Va1S16Nafmva1:2£J3S Val« 1 0Ncrm val :32- 52 
Valfll 7N0Jmval £2 39 Va1#17Ncrmva1:36.53 



VaW1SNe*m val 23 09 Val#18Ncrmval:40.54 
ValS1*Jonii V ali7.7 Val* 1 9N<srm val ;42,54 



Val«20Nc*m val' 29.3 7 Va.ff£C \ c — val -A2.5i 



Figure 16. (2) Eigenfunction correspondences 
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1-O.S. NoonVtl: .W.C372 t-0.73, NcnmVal: ,95.4»9 1-0. 75, NormW. ,54.3732 





MJ74.NormV a | 94.2924 1-0 72, NcrmVal- ,93,03i& 1-0,7, NonmVal; ,93,0124 





t=Oj63, NormVal: ,52.3462 UQJSS, NormVal: .5 1 .663 1 1*0.64. McrmV^: .51 .001 5 



# 



1-002 NormVal 90.3493 U0.0. NormVal: .4S.0333 



# 




Figure 17. (1 - 



— ) Family Tenth Eigenvalue Zoom 
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t-0.S,H«mV 4 l: ,56,71 07 




1-0 7&, NcrmVah ,55.9563 1-0.76, NwrmVil: ,55.217 





1=074, NcimVal: ,54*99 1-0 72, NciwiVah ,59.7964 1-0.7, NormVal ,53.1645 





t— 0.63, McrmVai: ,52.5966 t=0j66, NcrniVaJ: ,52.02)63 1*0.64, NcrmVal: ,51 .4663 



# 



1-6 .62. NormVal: 30:9151 1*0.6. HoimVal: 50.357$ 




FIGURE 18. (1 ) Family Eleventh Eigenvalue Zoom 
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55 

54.5 

1 54 

Sv53.5 
Lu 

53 



(a) 



SZ.5 



i t 



asa D.E9 

t VllUC 



□.71 Q.7Z 0.73 D.74 




0.735 0.73* 0.711 0.7+1 0.747 11.750 0.753 0.75S 0.759 
[ value 

(c) 





50.05 ■ 






«s 






50 - 


1 






49.95 - 




49.9 - 



0504 0.595 0.596 0.537 a.59S 0.599 0.6 0.ED1 0.50? 0.603 0.604 
e vain? 



Figure 19. Eigenvalue Dynamics Summary, (a) A collision. 

Tenth and eleventh (1 ) eigenvalues of H(t,-,-) as in Figs. 17 

and |18[ (b) A non-crossing tangency, with (analytically deter- 
mined) two-dimensional cigcnspacc. These are two (1 + +) eigen- 
values of the homotopy Go(t, ■, •) of Section[3]as in Figs. 36 and 37 



(c) A crossing, with (experimentally determined) two-dimensional 
eigenspace. These are the sixth and seventh (1-1 — ) eigenvalues of 
the homotopy Go(t, ■, •) of Section [3] as in Figs. 32 and 33 
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uoimmeo Eigenvalue Differences (a) 




e.?i5 o.nt hl.741 a:m 0.747 o.7» 0.753 0.73s n.7# 




-0.05 



□.534 D.SgS 0.595 0.597 0.59B 0.599 



0.60? 0.603 D.6W 



Figure 20. Eigenvalue Dynamics Summary, Eigenvalue Differ- 



ences Aj+i — Xj from Fig. 



19 



(a) A collision. 
) as in Figs 



Tenth and eleventh 
(b) A non- 



17 and 18 



(1 ) eigenvalues of H(t, 

crossing tangency, with (analytically determined) two-dimensional 
eigenspace. These are two (1 + +) eigenvalues of the homotopy 
Go(t, ■, •) of Section [3] as in Figs. 36 and 37 (c) A crossing, with 
(experimentally determined) two-dimensional eigenspace. These 

are the sixth and seventh (1-1 ) eigenvalues of the homotopy 

Go(t, ■, •) of Section [3] as in Figs. 32 and 33 
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Figure 21. G through G 3 , (1 + +) Eigenvalue Summary Plot 
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Figure 22. G through G 3 , (1 ) Eigenvalue Summary Plot 
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Figure 23. G through G 3 , (1 - +) Eigenvalue Summary Plot 
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Figure 24. G through G 3 , (1 + -) Eigenvalue Summary Plot 
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FIGURE 25. G through G 3 , (2) Eigenvalue Summary Plot 
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t-1. NflrmVal ,i7,£M& 

o 



M) 7, NoirnVal ,i 7.0639 




1.0.4, NarmVd: .16.6091 

a 




M) 5, N^mV*(; ,17,22* 1-0.3, N^mVtA ,17.1646 




1-0.3, NormVal: . 1 6.4 1 74 1^0 2, Norm Val: , 1 6.2£*6 




1-0, NonmVal-. ,16,0014 

m m 




Figure 26. G , (1 + +) EigVal 4 Homotopy 
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t-1. NwmVd ■ ,17.137* W.9,NwmVnl:,17.117 NwmVsl: .17.W6S 



0« 

«0 




WJ.7, N^mVtl: .17.065 1-0.6.. . N<*mV«A .1 7.074C t-0:5, ^ wm vd : ,17.1011 





1-0.4. <Jarm Val: . 1 7 . 1 4 1 4 1=0,3. Norm Ve* , 1 7. 1 E 1 8 1*0 ,£, Norm Va] : ,1 7.2 1 4 





*-0- 1 . N wm V* . 1 7 23S1 1-0, NonmVal , 1 7.2M& 





Figure 27. d, (1 + +) EigVal 4 Homotopy 
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t-1. Ne*mVd ,13,2054 W> 9, NwmVuf: ,15 5SW 1-0.3. NofmVuJ: ,13.S7M 



"1 S 11 




M) 7. NoimVal ,16.1674 1=0.6. NfxmVal: .16.4253 t-0-6, MormVsl: J 6.6425 



-I c 

1f 




1=0,4, NmmVal: , 1 6.32 1 7 1-0,3 . , NannVal : , 1 6J9609 1=0 2, Norm VbJ: , 1 7.0593 



1-0-1 . NnrmVal: .17.1 1 78 1-0, NonmVal , 17.1520 



Figure 28. G 2 , (1 + +) EigVal 4 Homotopy 
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U1 , mVaP: ,1 3 7958 M) .5, Ne™v«| ; .14 030i 1-D .3, MormVal : M £3i 5 




, NormVal:. ,15.110 1-0, NcrmVal- ,15.1*59 




^^^^ ^^^^ 

1-0.7, NwmVjft .144420 1-0.6, NsmiVtrf: .14,61 S4 t-0:6, Horn Val: ,14.7727 

W>.4, NmmVal: ,14.388* t-Q.3. MormVal: ,15.0(K 1-0 NormVal: ,15.0731 




Figure 29. G 3 , (1 + +) EigVal 4 Homotopy 
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1-1, NormVal ,17.£99& 1-0-9, NormVal: ,17.224 



313 

it 



A 



1*0 NormVal: ,17.1646 




Wfc5!t NormVal: , 1 7.06«9 1*0.6, NormVal : ,1 6.9363 





1-0 5. NormVal: ,16.7646 

w 



1*0.4, NormVal: ,16.8091 1*0.3, NormVal: ,16.4174 H) .2. NormVal: ,16.22,; 3 






U0. 1 , NormVal: , 1 6 .0668 1*0. NormVal: , 1 6.001 4 





FIGURE 30. G , (1 + +) EigVal 4 Nodal Homotopy 
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FIGURE 31. G , (1 + +) EigVal 4 Line Homotopy 
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Figure 32. (1 + -) EigVal 6 
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Figure 33. (1 + -) EigVal 7 
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Figure 34. (1 + +) EigVal 13 
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Figure 36. A = 1600 at t = .125, EigFcn 1 
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Figure 37. A = 1600 at t = .125, EigFcn 2 
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